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Direct capture in the ^^°Sn(n,7)^'^^Sn and ^^^Sn(n,7)^^^Sn reactions under r-process 

conditions 
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The cross sections of the ^^''Sn(n,7)^^^Sn and ^^■^Sn(n,7)^"'^Sn reactions are calculated in the direct 
capture model at low energies below 1.5 MeV. Using recent data from (d,p) transfer experiments 
on Sn and ^^^Sn, it is possible to avoid global input parameters with their inherent uncertainties 
and to determine all input to the direct capture model by local adjustments. The calculated direct 
capture cross sections of ^'^°Sn and "'^^'^Sn are almost identical and have uncertainties of less than 
a factor of two. The stellar reaction rates Na < av > show a slight increase with temperature. 
Finally an estimate for the influence of low-lying resonances to the stellar reaction rates is given. 
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Direct capture (DC) is expected to be the dominating 
reaction mechanism if the level density in the compound 
nucleus is low. This is typically found for light and/or 
neutron-rich nuclei, especially with magic proton or neu- 
tron numbers, at low energies which is the relevant energy 
range for nuclear astrophysics. Direct neutron capture 
has been identified experimentally for several stable tar- 
gets (e^. ^Li 0, [3, ^^O i, ^^O ^^Ne iQ, 

26Mg [iil, 4*^Ca [l3,[ll|), but it is obvious that neutron 
capture experiments are practically impossible for short- 
living radioactive targets like ^^^^Sn or I'^^gn^ Thus, the 
determination of the DC cross section for unstable tar- 
gets has to rely on theoretical predictions. 

The calculation of DC cross sections requires several 
ingredients. First of all, the electromagnetic transition 
must be well defined. This requires the transition energy 
Ej = E+Sn — Ex and thus the neutron separation energy 
Sn (or the masses of the target and residual nucleus) and 
the excitation energy E^ of the final state. In addition, 
spin and parity J'^ of the final state and its spectroscopic 
factor C^S are essential ingredients for the calculation. 
Finally, the DC cross section depends on the square of 
the overlap integral I 
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where O^^ is the electric dipole operator and u{r) and 
x(r) are the bound state wave function and scattering 
state wave function. These wave functions are calcu- 
lated from the two-body Schrodinger equation using a 
simple nuclear potential without imaginary part because 
the damping of the wave function in the entrance chan- 
nel by the tiny DC cross sections is very small [l^l- The 
present study is restricted to El transitions which are 
dominant in the DC cross section whereas higher multi- 
polarities like Ml or E2 are practically negligible if dom- 
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inant El transitions are allowed by the well-known elec- 
tromagnetic transition rules 13 
DC model can be found e.g. in 

The cross sections of the i'^°Sn(n,7)i'^iSn and 
i'^2Sn(n,7)i'^'^Sn reactions play an important role in r- 
process nucleosynthesis. In general, the influence of neu- 
tron capture cross sections on r-process nucleosynthesis 
is relatively small because under typical conditions an 
equilibrium between (n,7) and (7,n) reactions is reached. 
However, during freeze-out the cross sections become im- 
portant. This holds in particular for the i'^°Sn(n,7)-'^'^iSn 
reaction because of the larger neutron separation energy 
Sn = 5206 ± 13keV of ^-^^Sn compared to the smaller 
S„ = 2370±24keV for "3Sn (taken from the latest mass 
evaluation A detailed study of the r-p rocess nu- 

cleosynthesis around Aw 130 is given in [1^. The par- 
ticular importance of the ^■^°Sn(n,7)i"^iSn cross section 
is highlighted in [T6j . and the most important tempera- 
ture range is identified as 0.8 < Tg < 1.3 (where Tg is 
the typical notation for the temperature in 10^ K). This 
corresponds to thermal energies 70keV < kT < llOkeV. 
Because of the missing Coulomb barrier in neutron cap- 
ture, the stellar reaction rate per mol and unit volume 
TZ{T) = Na < crv > (the usual short term "reaction 
rate" will be used for TZ in the following) is mainly sensi- 
tive to the cross sections at energies around E « kT, and 
the temperature dependence of TZ{T) is small (for pure 
s-wave capture a 1/v and TZ(T) = const.). 

Up to now, the DC cross sections of the 
i^°Sn(n,7)i^iSn and i^2Sn(n,7)i^^Sn reactions have 
been calculated using global parametrizations of the 
required input parameters [l^, [ij [l^ • It was found that 
the DC cross section of the i^°Sn(n,7)i'^iSn reaction is 
very sensitive to the chosen parameters. At 30keV a 
variation over three orders of magnitude is found (see 
Fig . 9 of The recent (d,p) experiments on ^^^Sn 

and ^^^Sn allow for the first time to completely 
avoid global parametrizations. Instead, locally optimized 
parameters are used in this work for all ingredients of the 
DC calculation to minimize the resulting uncertainties. 
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I start with the analysis of the DC cross section for 
the doubly-magic ^^^Sn target nucleus. The bound state 
properties of the residual ^^'^Sn are well-known from the 
^^^Sn(d,p)^^'^Sn experiment and are summarized in 
Table H] The spectroscopic factors are compatible with 
unity (see Table I in [lo]); thus, 5 w 1.0 is adopted 
in the following calculations. Such large spectroscopic 
factors are expected for single-particle states above the 
doubly-magic ^^^^Sn. 



TABLE I: Properties of bound states in ^^^Sn and "^Sn (from 



19142311) and the considered El transitions. 



J" 


(keV) E (keV) 


C^S 


Vo (MeV) U 


Lf 








Sn 








3/2+ 





-5206 


0.10 


-39.40 


1,3^ 
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1/2+ 


332 


-4874 


0.10 


-40.04 
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5/2+ 


1655 


-3551 


0.10 


-36.82 


1,3^ 
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7/2" 


2628 


-2578 


0.70 


-47.30 


2,4^ 
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3/2" 


3404 


-1802 


0.70 


-46.97 


0,2^ 
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1/2- 


3986 


-1220 


1.00 


-45.66 


0,2-> 
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5/2- 


4655 


-551 


0.75 


-43.70 


2,4^ 
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133 


Sn 








7/2- 





-2370 


^ 1.0 


-46.51 


2,4^ 
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3/2- 


854 


-1516 


^ 1.0 


-45.93 


0,2^ 
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1/2- 


1363 


-1007 


^ 1.0 


-44.74 


0,2^ 
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5/2" 


2005 


-365 


^ 1.0 


-42.91 
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The nuclear potential is taken as the sum of a central 
and a spin-orbit potential 



Vir)^-Vofir)-VLs — fLS 

r ar 



with the central depth Vq, the spin-orbit strength Vls 
and the Woods-Saxon geometry 



f{r)= l+cxp(- 



(3) 



1 /3 

with the radius parameter R ^ Rq x Aj( , Rq = 1.25 fm, 
and a = 0.65 fm. 

In a first step the bound state wave functions u{r) are 
calculated by adjusting the depth Vq of the central po- 
tential (with Vls = 0) to the energy E <0 (see Table 
With an additional spin-orbit potential almost identical 
wave functions can be obtained using Vq = 45.5 MeV 
(45.0 MeV) and Vls = 18.6 MeV (22.0 MeV) for the 
bound L = 1 (L = 3) states. 

The second step is the calculation of the scattering 
wave function x(r). The optical potential can be ad- 
justed to experimental phase shifts for all partial waves 
or to the scattering length for the s-wave. Unfortunately, 
such data are not available for the unstable nuclei under 
study. As an alternative, the potential strength can be 
adjusted to the energies of single-particle states (as al- 
ready done for the bound states above). For light nuclei 
often a significant parity dependence for the potential 
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FIG. 1: (Color online) DC cross section of the 

^''^Sn(n,7)^'^'*Sn reaction. The contributions of the bound 
states in Table [l] are shown by colored dashed and dotted 
lines. The full black line represents the sum over all bound 
states. 



depth Vq is found. However, with increasing mass this de- 
pendence decreases, and e.g. for ^^Ca (above the doubly- 
(2) magic ^*Ca) it is found that Vq of the bound L = 1 states 
(derived from the bound state energies) and Vq of the s- 
wavc (derived from the scattering length) agree within 
about 1 %. (This result is independent of details of the 
geometry of the potential; also for a folding potential the 
deviation is only about 1%.) [T^j. Because of the mi- 
nor difference of Vq and Vls for the L = 1 and L — 3 
bound states, I adopt the average of Vq = 45.3 MeV and 
Vls = 20.3 MeV for the calculation of the scattering wave 
functions x('')- 

Now all parameters for the calculation of the overlap 
integrals I in Eq. ^ are fixed by local adjustments to 
properties of ^'^^Sn = ^'^^Sn (g) n, and the DC cross sec- 
tions can be calculated without any further adjustments 
or parameters from global studies. The result for the 
^^^Sn(n,7)^'^'^Sn cross section is shown in Fig. [T] A dis- 
cussion of uncertainties will be given later. 

Exactly the same procedure is repeated for the 
^•^°Sn(n,7)^'^^Sn reaction. The bound state properties 
of the L = 1 and L = 3 bound states are taken from 
the recent ^■^°Sn(d,p)^'^^Sn experiment [l^. Very similar 
to ^'^•^Sn, no fragmentation of the levels has been found 
in ^■^^Sn which is somewhat unexpected for the semi- 



magic '^^Sn core (compared to the doubly-magic Sn 
core in the previous case). The resulting average param- 
eters Vo = 46.2 MeV and Vls = 21.1 MeV are derived 
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FIG. 2: (Color online) DC cross section of the 

^^''Sn(n,7)'^'^'^Sn reaction. The contributions of the bound 
states in Table [l] are shown by colored dashed and dotted 
lines. The full black line represents the sum over all bound 
states. 



from Vb = 46.6 MeV (45.8 MeV) and Vls = 20.3 MeV 
(21.9 MeV) for the i 1 (L = 3) bound states. The 
potential parameters remain very close to the data for 
^^■^Sn and confirm the similarity of ^'^"'^Sn and ^'^'^Sn. 

The bound states with even parity are characterized 
by a particle-hole structure [l^. Thus, they have much 
smaller spectroscopic factors. These states are not suited 
for a determination of the potential depth Vq which shows 
a broader spread. A spectroscopic factor of C^S = 0.1 
has been assumed for these states which is in agreement 
with the upper limit of w 0.3 given in [l^ but somewhat 
lower than the average value of 0.347 for compiled spec- 
troscopic factors pJi, il^l- The DC cross sections for the 
bound states with even parity are much smaller than for 
the odd-parity bound states. The total DC cross section 
(summed over all transitions) does not depend strongly 
on the assumed value of C^S = 0.1 for the weak transi- 
tions to the bound states with positive parity (see Fig. [2]). 

From the DC cross sections in Figs. [T] and [2] stellar re- 
action rates TZ{T) = Na < cfv > can be calculated. Note 
that the laboratory reaction rate Ti\ah and the stellar re- 
action rate TZ* are practically identical in the important 
temperature range around Tg « 1 [2^. The reaction 
rates of both reactions under study are very similar and 
show a weak temperature dependence (see Fig. The 
results can be simply parametrized by a three-parameter 
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FIG. 3: (Color online) Stellar reaction rate 7^(T) = 
Na < (JV > for the i^°Sn(n,7)i^iSn (red dotted) and 
^^^Sn(n,7)^^^Sn (blue dashed) reactions. 



parabolic fit 

7^(T) = NA<(yv >? 



(ao + aiTg + a2Tg ^ 



s mol 



(4) 



with ao = 16811 (16321), ai = 2291 (2236), and aa = 
700 (870) for "Ogn (i32Sn). The deviations of the fit are 
1 — 2 % over the full temperature range under study. 

Uncertainties of the DC cross sections are studied by 
a variation of the different parameters of the calculation 
within reasonably estimated ranges and by considering 
the experimental uncertainties of the bound state prop- 
erties. The uncertainty of the neutron separation ener- 
gies Sn and the excitation energies lead typically to 
uncertainties for the transition energy of less than 
10%. Together with the E'^ dependence of the El tran- 
sition strength a typical uncertainty of about 10 — 30 % 
is found for the various transitions under study. A vari- 
ation of the potential geometry (using a larger value of 
i?o = 1.4fm instead of i?o = 1.25fm) and readjusting 
the potential depths leads to variations of the DC cross 
section between 10 — 20%. A reduction of the potential 
depth Vq by 3% reduces the DC cross section by about 
15%. The spectroscopic factors C'^S have uncertainties 
of about 30% which enter linearly into the DC calcula- 
tion. Combining all the above uncertainties of the order 
of 10 — 30%, a total uncertainty below 50% is a rea- 
sonable estimate for the total DC cross section of the 
"0Sn(n,7)"iSn and i32Sn(n,7)"3Sn reactions. 

For the ^■^^Sn(n,7)-'^'^'^Sn reaction reaonable agreement 
with the three predictions in [l7| is found whereas the 
new result is lower by a factor of slightly above 2 (slightly 
below 2) than the calculation in |3 ([isj). The energy 
dependence of all calculations jl3l |l8| is very similar be- 
cause it is essentially defined by the angular momenta 
in the entrance channel in combination with the electro- 
magnetic selection rules. 

The obtained results for the ^■^°Sn(n,7)^'^^Sn reaction 
are slightly below but very close to the calculations shown 
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in Fig. 4 of [l9[. This is not surprising because the same 
bound state properties {J'" and Er^) are used. The essen- 
tial difference between this work and 1191 is the replace- 
ment of the global optical potential in [1^ by the locally 
optimized potential which reduces the uncertainties for 
the calculated aoc- 

The new aoc for ^^°Sn(n,7)^^^Sn is about a factor 
of two below the highest result by Rauscher et al. [l7j . 
There are two further calculations in 171 with much 



smaller cross sections which result from the fact that 
some of the bound states in Table U are unbound in the 
corresponding calculations. However, the dramatic re- 
duction of the DC cross section in [T^] is an artifact from 
the separate treatment of the entrance and exit chan- 
nels. If the L = \ bound states were indeed unbound, 
the L = \ strength would be located close above thresh- 
old and show up as resonances in ajjc (and increase ctdc 
via transitions to bound positive-parity states in ^'^^Sn 
instead of reducing aoc)- This can be simulated by a 
reduction of the potential depth Vb, but is not taken 
into account in [17| using a fixed potential in the en- 
trance channel. E.g., using Vq = 41.0 McV (instead of 
46.2 McV) leads to a strong 3/2~ resonance at about 
73keV with a total width F w 58keV and a total cross 
section of 4.2 mb in the resonance maximum, i.e. a factor 
of about 50 higher than the standard calculation shown 
in Fig. [5J The resulting stellar reaction rate TZ becomes 
temperature-dependent and would be a factor of 10 — 20 
higher than the result in Fig. [3] because of this artificial 
3/2^ resonance. However, such a strong resonance has 
been excluded by the transfer data [l9| . 

Finally, predictions of the ^^°Sn(n,7)^^^Sn and 
^'^^Sn(n,7)-^'^'^Sn cross sections from the statistical model 
have to be discussed briefly. As pointed out e.g. in psj . 
the statistical model is not applicable below Tg « 1.4 
for "2sn and below Tg « 0.2 for ^^°Sn because the 
level density is too low. The limit for ^"^^Sn may even 



be higher if one takes into account that surprisingly low 
fragmentation of strength and very similar properties of 
^■^^Sn and ^^'^Sn were found in the transfer experiments 
p^ . [20j . As a consequence, large deviations are found 
for predictions from the statistical model using different 
ingredients (for details sec Fig. 1 of [l5| and discussion). 
Thus, a better estimate for resonant contributions might 
be the procedure of lowering the potential depth Vb (as 
outlined above). From the spectroscopic factors in [1^ 
(see also Table |T]) the missing w 25 % of the L = 1 or 
L = 3 strengths may be located above threshold, but 
below the detection limit of [l^. A resonance with full 
3/2~ strength would lead to an enhancement of the stel- 
lar reaction rate 7?. by a factor of 10 — 20; thus, a weaker 
resonance with 25% of the strength should enhance TZ 
not more than a factor of 2.5 — 5 if located close above 
the threshold, and the resonant enhancement is decreas- 
ing for higher-lying resonances. Such an enhancement is 
only expected in the ^■^°Sn(n,7)^'^^Sn reaction, but not for 
the ^'^^Sn(n,7)^'^'^Sn reaction because there are no bound 
states with positive parity in i33Sn mi,!!!. 

In summary, the direct capture cross section of the 
^^°Sn(n,7)^^^Sn and ^^^Sn(n,7)^^^Sn reactions has been 
calculated using local parameters which could be derived 
mainly from recent (d,p) transfer experiments [l^, [l^] ■ 
The DC cross sections of -'^'^^Sn and ^"^^Sn are almost 
identical and could be determined with relatively small 
uncertainties of less than a factor of two. Additional reso- 
nant contributions may enhance the stellar reaction rate 
by up to a factor of 5 for ^'^'^Sn depending on whether 
the remaining L = 1 and L ~ 3 strength is located in 
a narrow energy window close above threshold. Huge 
enhancements of the reaction rate 7?. of a factor of 10 
or even 100 (as discussed in [3) are excluded by the 
present study. 

This work was supported by OTKA (NN83261). 



[1] Y. Nagai, M. Igashira, T. Takaoka, T. Kikuchi, T. Shima, 
A. Tomyo, A. Mengoni, T. Otsuka, Phys. Rev. C 71, 
055803 (2005). 

[2] T. Kikuchi, Y. Nagai, T. S. Suzuki, T. Shima, T. Kii, M. 
Igashira, A. Mengoni, T. Otsuka, Phys. Rev. C 57, 2724 
(1998). 

[3] M. Igashira, Y. Nagai, K. Masuda, T. Ohsaki, H. Ki- 

tazawa, Astrophys. J. 441, 89 (1995). 
[4] J. Meissner, H. Schatz, J. Gorres, H. Herndl, M. Wi- 

escher, H. Beer, F. Kappeler, Phys. Rev. C 53, 459 

(1996). 

[5] T. Ohsaki, M. Igashira, Y. Nagai, M. Segawa, K. Muto, 

Phys. Rev. C 77, 051303(R) (2008). 
[6] H. Beer, P. V. Sedyshev, W. Rochow, P. Mohr, H. Ober- 

hummer, Nucl. Phys. A705, 239 (2002). 
[7] A. Tomyo, Y. Nagai, Y. Nobuhara, T. Shima, H. Makii, 

K. Mishima, M. Igashira, Nucl. Phys. A718, 527 (2003). 
[8] P. Mohr, H. Beer, H. Oberhummer, G. Staudt, Phys. 

Rev. C 58, 932 (1998). 



[9] P. Mohr, H. Beer, H. Oberhummer, W. Rochow, 
P. V. Sedyshev, S. Volz, A. Zilges, Phys. Rev. C 60, 
017603 (1999). 

[10] H. Beer, C. Coceva, P. V. Sedyshev, Yu. P. Popov, 

H. Herndl, R. Hofinger, P. Mohr, H. Oberhummer, Phys. 

Rev. C 54, 2014 (1996). 
[11] P. Mohr, H. Oberhummer, H. Beer, W. Rochow, V. KoUe, 

G. Staudt, P. V. Sedyshev, Yu. P. Popov, Phys. Rev. C 

56, 1154 (1997). 
[12] E. Krausmann, W. Balogh, H. Oberhummer, T. 

Rauscher, K.-L. Kratz, W. Ziegert, Phys. Rev. C 53, 

469 (1996). 

[13] Yi Xu and S. Goriely, Phys. Rev. C 86, 045801 (2012). 
[14] http://www-nds.iaea.org/amdc/, version amdc.114; G. 

Audi and W. Meng, private communication. 
[15] R. Surman, J. Beun, G. C. McLaughlin, W. R. Hix, Phys. 

Rev. C 79, 045809 (2009). 
[16] J. Beun, J. C. Blackmon, W. R. Hix, G. C. McLaughhn, 

M. S. Smith, R. Surman, J. Phys. G 36, 025201 (2009). 



5 



[17] T. Rauscher, R. Bieber, H. Oberhummer, K.-L. Kratz, 

J. Dobaczewski, P. MoUer, M. M. Sharma, Pliys. Rev. C 

57, 2031 (1998). 
[18] S. Chiba, H. Koura, T. Hayakawa, T. Maruyama, T. 

Kawano, T. Kajino, Phys. Rev. C 77, 015809 (2008). 
[19] R. L. Kozub et ai, Phys. Rev. Lett. 109, 172501 (2012). 
[20] K. L. Jones et ai, Phys. Rev. C 84, 034601 (2011). 
[21] Online database ENSDF, 

http://w'ww. nndc.bnl.gov/ensdf/. 



[22] Y. Khazov, I. Mitropolsky, A. Rodionov, Nucl. Data 

Sheets 107, 2715 (2006). 
[23] Y. Khazov, A. Rodionov, F. G. Kondev, Nucl. Data 

Sheets 112, 855 (2011). 
[24] S. Goriely, Phys. Lett. B. 436, 10 (1998). 
[25] T. Rauscher and F.-K. Thielemann, At. Data Nucl. Data 

Tables 75, 1 (2000); http://www.nucastro.org. 



